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ABSTRACT 

We present results of our calculations of NLTE model stellar atmospheres for hot Population III stars composed of 
hydrogen and helium. We use our own computer code for calculation of spherically symmetric NLTE model atmospheres 
in hydrostatic and radiative equilibrium. The model atmospheres arc then used for calculation of emergent fluxes. 
These fluxes serve for the evaluation of the flow of high-energy photons for energies higher than ionization energies of 
hydrogen and helium, the so-called ionizing photon fluxes. We also present the time evolution of the ionizing photon 
fluxes. 

Subject headings: Stars: atmospheres - radiative transfer - stars: Population III - first stars 



1. INTRODUCTION 

At the beginning of the era of modern NLTE stellar at- 
mosphere modelling at the early 70's, the chemical com- 
position of computed model atmospheres had to be as 
simple as possible to enable their calculation at comput- 
ers that were available at that time. Therefore, these 
models included only minimum number of elements, se- 
lected according to their importance for the atmospheric 
stru cture properties. Firs t mod els were purely hydro- 
gen (lAuer fc M ihalas 1969, 19701), later also helium was 
~ Mihalas fe AuerT 119701 1 Auer fc Mihalal H97l 
19 7311 19 761) to calculations. Due to computer 
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limitations, the absorption by heavier elements was not 
considered, and, of course, also the absorption by a 
huge number of metal lines in the ultraviolet spectral 
region was omitted. These metal lines cause lowering 
the ultraviolet flux and enhancing the flux in the vi- 
sual region. This important effe ct is known as line blan- 
ket ing (see, e.g., iMihalasI 119781) . Great effort (started 
by lAndersoiil 119851 119891) was necessary to include the 
line blanketing into NLTE model atmospheres calcula- 
tions, and nowadays it is possible to calculate NLTE 
model atmospheres with this effect taken into account 
properly. One of the most widely used codes for cal- 
culation of static NLTE line blanketed plane-parallel 
model atmosphe res is the code TLUSTY ()Hubenvl H98H; 
iHubenv fc Lanzlll995[ ). NLTE line-blanketing can be also 
treated by other model at mosphere codes, like the model 
atmosphere package TMAP ([Werner et all 120031) . For the 
case of NLTE mo ving atmospheres wit h line blanket- 
ing, codes CMFGEN (iHillier fc Millerlll998D. PHOENIX (e.g. 



Hauschildt et al 
1997t IPuls et al 



1997), FASTWIND tiSantolava-Rev et all 
2005), and WMBasic (e.g. IPauldracbl 



2003L and references therein) are available. 



All these mentioned model atmosphere codes were de- 
veloped and used with the primary objective to repro- 
duce the observed emergent radiation from hot stars. 
The hyd rogen-helium sta tic model atmospheres calcu- 
lated by lAuer fc Mihalai served only as a first approxi- 
mation of NLTE hot stellar atmosphere modelling. How- 
ever, there is a group of stars, whose chemical com- 
position corresp onds exactly to what was calculated by 
lAuer fc Mihalasl in the early 70's, namely the first stars 
in the Universe. Since the Universe at its beginning was 



practically free of metals and consisted only of hydrogen, 
helium, a nd a very little amount of lithium (see, e.g. , 
reviews of IBromm fc Larson! [200l I Johnson et al. 2008; 
iBromm et al.ll2009D . the first stars may be safely assumed 
to consist only of hydrogen and helium, since the influ- 
ence of lithium on the stellar atmospheric structure is 
negligible. 

Thus, in order to obtain model atmospheres of first 
stars, it is simply sufficient to use the first codes for 
calculation of NLTE model atmospheres from the early 
70's, or to repeat such calculations with contemporary 
sophisticated model atmosphere codes and just to skip 
metals from the ca l culati ons. The latter was done by 
iTumlinson fc Shulll (|2000l) . who used the code TLUSTY 
to calculate simple hydrogen- helium model atmospheres 
and rediscovered that "A model atmosphere is necessary 
because a simple blackbody curve for each T e g will not 
accurately reproduce the spectrum ... " (cf. Unsoldi ri955l . 
where the history of the finding that the stellar radiation 
is not the blackbody one is thoroughly descibed). They 
concluded that neglecting all the metal opacities leads to 
enhancement of the UV flux emerging from the model 
stellar atmosphere, a fact that was well known before, 
but in opposite formulation, namely that adding metal- 
lic opacity reduces the ultraviolet stellar flux. 

The presence of high UV flux emerging from the 
hydrogen-helium model atmospheres describes an im- 
portant source of radiation, which is able to ionize the 
intergalactic medium in the early Universe evolution 
( Venkates an et alJ 120031 ) . Indeed, high fluxes at high 
frequencies caused by missing stellar opacity in those 
spectral regions are extrem ely important in cosmology 
(|Venkatesan fc Trur an 2003) . Model atmospheres of first 
stars and their infl uence on ionizing radi ation were also 
studied in detail by iSchaererl (|2002l 120031 ) . 

The absence of metals in the atmospheres of first 
stars has an additional effect. Since the winds of hot 
stars are driven by radiation scattered by spectral lines, 
there was a question whether the first stars have stel- 
lar winds. Existence or non-existence of the stellar 
wind has dramatic influence on the fi rst star evolu- 
tion. IBromm. Kudritzki. fc Loebl (|2001l ) calculated uni- 
fied model at mospheres (using an improv ed version of 
the method of iSantolava-Rev et al.l fl997t ) and did not 
find any effect of the velocity field on line profiles for low 
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mass-loss rates. iKudritzkH (|2002f) tested the existence of 
winds of first stars by calculating wind models using the 
depth dependent line force multiplier parameters fc, a, 
and 5. He tried to mimic the zero metallicity of the first 
stars using a method, which implicitly assumes driving 
of the wind by metalic lines. He pointed out the possi- 
bility of multicomponent effects in t he wind of first stars . 
Thes e were later studied in detail bv lKrticka et al.l (2003, 
2010). However, the principal question of existence or 
non-existence of the winds of first stars deserve d more 
detailed study. The code of IKrticka fc Kubatl ()2004f ) 
enabled calculation of the radiative force directly using 
actual level populations and opacities of individual ele- 
ments at each depth point, without the necessity to in- 
troduce parameterization of the radiative force by means 
of the line force multipliers. Consequently, this method 
gives implicitly depth-dependent line force and allows to 
include co nsistently any chemical c ompo sition. Using 
this code, IKrticka fc Kubatl (|2006l 120091) studied stel- 
lar winds of firs t stars and confirmed the fact found by 
iKudritzkl (|2002f ) that the first stars have extremely weak 
winds. These calculations also showed that using static 
approximation for calculation of model atmospheres of 
first stars is adequate. For these calculations, which were 
done using the core-halo approximation, it was necessary 
to determine the flux at the lower boundary of the wind 
region. To this end, the emergent radiation from the 
static spherically symm etric NLTE model atmospheres 
was used. Consequently, Krticka fc Kubatl did not check 
if the plane-parallel approximation, which has been fre- 
quently used by others, is adequate for modelling the at- 
mospheres of fir st stars. Since th e NLTE atmospheric 
models used by IKrticka fc Kubatl have not been pub- 
lished anywhere in detail, in this paper we present an 
extended set of NLTE model atmospheres used for this 
purpose. 

2. MODEL ATMOSPHERES 
2.1. Validity of static approximation 

For the case when the thickness of the atmosphere is 
small compared to the stellar radius, the plane-parallel 
approximation offers a good choice. This happens for 
dwarf stars with relatively high surface gravitational ac- 
celeration (logg > 4). However, small deviations from 
the plane- parallel approximation were fo und even for 
subdwarfs (jGruschinske fc Kudri tzki 1979J). For the case 
of white dwarfs they were found only in cores of som e 
strong lines and were usually very small (Kubat 1995). 

On the other hand, for Population I and II stars with 
lower surface gravities, giants and supergiants, the at- 
mosphere is more extended, and it can no longer be con- 
sidered as thin with respect to the stellar radius. In 
such case a spherically symmetric approximation is more 
realistic than the plane-parallel one. However, due to 
a weaker gravitational force in the atmospheres of gi- 
ants and supergiants, the radiation force may easily over- 
come gravity there and give rise to a stellar wind. Con- 
sequently, such atmosphere can not be considered as a 
static one. Winds in hot star atmospheres are driven by 
radiation force generated by momentum transfer from 
absorptions and scatterings in lines of heavy elements, 
with a help of the force generated by continuum transi- 
tions and electron scattering. Note that the force gen- 



erated by continuum transitions and electron scattering 
alone is not sufficient to drive the wind. 

The first stars in Universe were born from pure 
hydrogen- helium mixture with no heavier elements. Con- 
sequently, since the radiation force does not overcome 
grav ity there, the line driven wind does not exist there 
(see IKrticka fc Kubatl 120061 I2009D and atmospheres of 
these stars may be described using static approximation 
in both plane-parallel and spherically symmetric geome- 
try. 

2.2. Method of calculation 

Our method for calculation of NLTE model stellar at- 
mospheres is based on an accelerated lambda iteration 
( ALI) method and has been descr ibed in some detail 
in IKubatl (I19941 . 119961 Il997bl . 120011 ) and summarized in 
IKubatl (|2003f l . Using our code we calculated a grid of 
hydrogen-helium spherically symmetric NLTE model at- 
mospheres in hydrostatic and radiative equilibrium. The 
helium abundance Yne = ^He/«H = 0.1. We consid- 
ered a 16 level hydrogen atom (15 levels of Hi + 1 level 
of H n) and a 50 level helium model atom (29 levels of 
He I + 20 levels of Hen + 1 level of Hem). Individual 
hydrogen levels correspond to main quantum numbers 
n = 1,...,15. For He I, all levels up to n — 4 were 
taken separately according to orbital quantum number 
I, for 5 < n < 9 we took 2 averaged levels for each n, 
one for singlets, one for triplets. The levels of Hen were 
considered similarly as the hydrogen levels, however for 
n = 1, . . . , 20. Details of the model atoms used (oscillator 
strengths, photoionization cross sections, collisional ex- 
ci tation and ionizatio n rates, line profiles) can be found 
in IKubat et all (fl999h . 

For temperature structure determination we used a 
combination of three methods. At large depths below 
the stellar radius i?* = r(rn w 2/3) (tr is the Rosseland 
optical depth) it was the differential form of the radia- 
tive equilibrium equation, above the stellar radius the 
integral form of radiative equilibrium equation, and for 
outermost layers the el ectron thermal balance method 
was u sed (for details see IKubatl IT9961 . 120011 : IKubat et all 
fl999t ). 

3. RESULTS OF CALCULATIONS 
3.1. Model atmosphere parameters 

Basic parameters (luminosity L*, mass M*, and ra- 
dius i?*) of the first stars are taken from evolution- 
ary calculations. For t he original calculations used in 
IKrticka fc Kubatl (120091) a limited se t of model param- 
eters (taken from lMarigo et all 120011 ) was used. In ad- 
dition to them, here we extended the set by using addi- 
tional parameters from M arigo et all and we also used 
other sources of stellar parameters following from other 
evolutionary calculations. Values of masses of these 
stars are still matter of debate and they are different 
in different calculations (cf. d ifferent results obtain ed 
by evolutionary calcu lations of lYqon fc Langerl 120051 or 
lYoon et al.ll2006l and lHirschill2007f ). 

Firs t stars model atmospheres calculated by iSchaererl 
(2002) cover stellar masses from 5 M Q up to extremely 
high masses of 1000 M Q . He used the code TLUSTY for 
calculation of H-He model atmospheres for stars with 
Toff > 20000 K, for colder stars he adopted the Ku- 
rucz LTE model atmospheres with a very low metallicity. 
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TABLE 1 

Model atmosphere parameters and and stellar ionizing photon fluxes Qj for s elected model atmospheres of first stars 

WITH PARAMETERS FROM SCHAERER (2002). 
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200 400 600 800 1000 1200 1400 
MA) 

Fig. 1. — Plot of ultraviolet fluxes fro m selected mod el atmo- 
spheres of first stars with parameters from Schacrcr (2002). Labels 
correspond to the first column of the Table [T] 

We calculated spherically symmetric NLTE mo del atmo- 
sphere s for parameters given in the Table 3 of iSchaererl 
(2002), these parameters are listed here in the Table [TJ 
corresponding temperature structures as a functio n of 
the column mass depth m are shown in the Figure IA11 
The column mass depth is related to the Rosseland op- 
tical depth by a relation 

dm= (i?2/r 2 )(p/ XR )dr R , (1) 

where p is the density and xb. t ne Rosseland opacity. All 
models display a rise of temperature in the outer parts of 
the atmosphere, where the strongest lines are optically 
thick and the continuum is optically thin. This is a typ- 
ical NLT E effect in model atm ospheres without metals 
(see, e.g. lAuer fc Mihalaill970D . In model atmospheres 
with metal line b lanketing included, thi s temperature rise 
disappears (e.g. JLanz &: Hubenyll2007t ). 

To check the validity of the plane-parallel approxima- 
tion for the first stars by direct comparison of spherically 
symmetric and plane-parallel models calculated with the 
same code (which has not been done yet), we also calcu- 
lated NLTE plane-parallel model atmospheres for these 
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MA) 

Fig. 2. — Plot of ultraviolet fluxes from selected model atmo- 
sphere s of zero-age main sequence first stars from IMarieo et al.l 
(2001). Labels correspond to the first column of the Table [2] 

stellar parameters. We found that the difference between 
the temperature structures of spherically symmetric and 
plane-parallel model atmospheres is relatively small, sim- 
ilarly to other ca ses where the atmospheric extension is 
not too larg e (cf. | Gruschinske fc Kud ritzki 1979t iKubatl 
119961 H997ai l3. TT999T ). Also the difference between emer- 
gent fluxes from spherically symmetric and plane-parallel 
model atmospheres (the sphericity effect) is small, the 
flux from the spherically symmetric model atmosphere is 
slightly lower. 

Although it was initially assumed that Population III 
stars shoul d have very la r ge ma sses of the order of hun- 
dreds M Q , IMarieo et all |2001f ) calculated an extensive 
grid of evolutionary models of Population III stars with 
masses 0.7 M© < M < 100 M Q with the assumption of 
zero mass-loss. 

For selected masses from this grid we calculated NLTE 
model atmospheres for models of stars with zero age. 
Their param eters (L*, M*, R*) w ere taken from the on- 
line Table in iMarigo et al.l ()2001[ ) and are listed in the 
Table [2j Temperature stru ctures for selected models are 
plotted in the Figure IA2I All models display the same 
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TABLE 2 

Model atmosphere parameters and and stellar ionizi ng photon fluxes Qj for selected model atmospheres for parameters 

FROM IMARIGO ET AL.I (120011) . 
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typ e of the rise of temperature as the models in the Fig- 
ure EU 

3.2. Ionizing fluxes and their evolution 

As an immediate by-product of the model atmosphere 
calculation we can obtain the amount of radiation es- 
caping from the stellar photosphere. This radiation sig- 
nificantly influences the circumstellar matter and it may 
be the reason for stellar winds. Fluxes from m ode l at- 
mospheres listed in the Table [2] (see also Fig. IA2I) for 
selecte d parameters from the large set of iMarigo et al.l 
(2001) were use d as a lower boun d ary co ndition in wind 
calculations by iKrticka fc Kubat] (|2009f) for analysis of 
the existence of stellar winds. 

For the case of Population III stars, the amount of 
highly energetic emergent radiation is of extreme impor- 
tance, because it influences the ionization state of in- 
terstellar and intergalactic matter (therefore this radia- 
tion is often referred to as the ionizing radiation) and, 
consequently, it is importan t for the cosmic evolution 
IjVenkatesan & Truranll20M ). 

Ultraviolet emergent fluxes from model atmospheres 
selected from the lists in the Tables [T] and [5] are plotted 
in the Figures [T] and [H respectively. It is evident that 
the hottest and most massive (and hence most luminous) 
stars provide maximum flux in the far-ultraviolet region, 
which is the extremely important spectral region for gen- 
eration of the "ionizing flux" . Contribution to this flux 
from the low-massive zero-age main sequence first stars 
is several orders of magnitude weaker. 

For all models we calculate the stellar ionising photon 
flux Q (in photons s" 1 ) as 



Q = 



F » A 
— dis 

hv 



(2) 



(z^ion is the ionization frequency of the particular ioniza- 
tion edge) for Hi, Hei, and Hen ionization edges. The 
values of the stellar ionizing fluxes are listed in the Ta- 
bles [1] a nd [2] Comparin g our results with those calcu- 
lated by iSchaererl (|2002l ) by a totally independent code 
(TLUSTY) we find a g ood agree ment. This supports the 
conclusion that both iSchaereil 's and our codes give cor- 
rect results. 

We also plot the dependence of the number of ionizing 
photons for H I, He I, and He II ionization edges on stellar 
luminosity in the Figure [3l As expected, the number of 
ionizing photons rises with the stellar luminosity, the rise 
is significantly steeper for L < 3.5Lq. 

In addition, we studied the evolution of ionizing fluxes 
with the stellar evolution. For stellar masses 10 M Q , 
20 M Q , 50 M q and 1000 M Q we calculated model atmo- 
spheres for evolved Population III stars. Parameters of 
these models are listed in Tables iBll - IB4I 

The ultraviolet emergent fluxes of some of the calcu- 
lated evolved model atmospheres for masses 100 M Q and 
50 M are shown in the Figured] The latter figure shows 
that the high ionizing flux present at zero-age main se- 
quence first stars decreases with their evolution, as the 
stellar effective temperature decreases. 

In order to describe the changes of the ionizing flux 
with time qua ntita t ivel y, we calculated for each model 
listed in Tables [BT1 - |B4l the ionizing photon flux after the 
Eq. ^ for each ionization (Hi, Hei, Hell), and we listed 
these numbers in corresponding tables. The changes of 
these fluxes with the stellar evolution are shown in Fig- 
ure The fluxes first rise within the magnitude order of 
the initial value, then after reaching the maximum value 
they quickly drop by several orders of magnitude. The 
maximum value flux depends on the stellar mass, it is 
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Fig. 3. — Number of ionizing photons Q as a function of the 
stellar luminosity for H I (upper panel) , He I (middle pa nel), and 
Hell ( lower p anel) continua. Stella r parameters are after lSchaererl 
1(2003 . +) and lMarigo et al.l l(200ll . 0). 



reached sooner for more massive stars. The relative im- 
portance of first stars for the generation of ionizing pho- 
tons at different ages is shown in the Fig. [5] The massive 
stars contibrute by several orders of magnitude more, 
but since they evolve faster, their contribution lasts for 
much shorter time period than for stars with lower mass. 
For example, the 100 M star stops its contributing after 
3 • 10 6 years, while the 10 M Q star contributes for about 
19 • 10 6 years. However, the total contribution of the 
10 M Q is smaller. 
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Fig. 4. — Plot of ionizing fluxes from selected model atmospheres 
after IMarigo et"aLl ||2001|T for 100M Q (upper panel) and 50M Q 
(lower pa nel) f or d ifferent ages of first stars with parameters given 
in Tables lBll - |B4l . The curves are labeled with the stellar age in 
megayears. 

4. CONCLUSIONS 

We presented results of the NLTE model atmosphere 
calculations of first stars (Population III stars) assuming 
hydrogen-helium composition. These model atmospheres 
represent an extended set of model atmospheres used as 
a lower boundary flux condition i n win d analyses of first 
stars bv lKrticka ,& Kubat] (|2006l I2OO90 . For the models 
we also calculated emergent fluxes and, in addition, the 
"ionizing photon fluxes", which collect the contribution 
of the radiation in the ultraviolet part of the spectrum. 
Photon fluxes were calculated for radiation with shorter 
wavelength than Hi, Hei, and Hen ionization edges. We 
also studied the time evolution of all three photon fluxes 
and we found strong dependence on the stellar mass, the 
more massive stars are more important contributors to 
ionizing photon flux. This underlines the importance of 
the question of the maximum stellar mass for zero age 
first stars. 
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A. TEMPERATURE STRUCTURE OF CALCULATED MODEL ATMOSPHERES 
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Fig. Al. — Run of temperature in NLTE model atmospheres of first stars with parameters after Schaerer (2002) listed in Table [T] (labels 
correspond to the first column of the table). The independent variable m is the column mass-depth (see Eq. [TJ. 



NLTE model atmospheres of first stars 



9 



100000 



90000 



80000 



70000 



60000 



50000 



40000 



30000 



20000 



10000 




1e-06 0.0001 0.01 



100 10000 



log m (g cm ) 



Fig. A2. — Run of temperature in NLTE model atmospheres of ZAMS first stars with stellar parameters after Marieo ct al. (2001). 
Individual labels Mxxx correspond to to the first column of the Table |21 The independent variable m is the column mass-depth (sec Eq.[Tj. 
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TABLE Bl 

Model atmosphere parameters and and stellar ionizing photon fluxes Qi f or selected model at mospheres from the 

EVOLUTIONARY SEQUENCE FOR IOOMq FOR PARAMETERS FROm IMARIGO ET AL.I (120011 ). 
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2 


91 


10 47 


2 


874 


10 6 


30.3 


100 


2 


218 


10 6 


40460 


3 


475 


1 


00 


10 49 


2.47 


10 48 


2 


77 


10 47 


2 


889 


10 6 


37.4 


100 


2 


239 


10 6 


36470 


3 


291 


9 


09 


1Q 48 


1.74 


10 48 


1 


55 


10 47 


2 


900 


10 6 


42.5 


100 


2 


249 


10 6 


34280 


3 


181 


8 


40 


1Q 48 


1.35 


10 48 


1 


05 


10 47 


2 


911 


10 6 


48.0 


100 


2 


259 


10 6 


32280 


3 


075 


7 


67 


10 48 


1.01 


10 48 


6 


74 


10 46 


2 


920 


10 6 


53.4 


100 


2 


265 


10 6 


30620 


2 


982 


6 


93 


10 48 


7.31 


10 47 


4 


25 


10 46 


2 


925 


10 6 


56.4 


100 


2 


280 


10 6 


29850 


2 


935 


6 


58 


10 48 


6.18 


10 47 


3 


35 


10 46 


2 


975 


10 6 


125. 


100 


2 


311 


10 6 


20090 


2 


241 


5 


86 


10 47 


6.01 


10 42 


3 


15 


1Q 39 


3 


014 


10 6 


510. 


100 


2 


327 


10 6 


9977 


1 


022 


3 


86 


10 42 


1.05 


10 36 


1 


32 


10 27 



B. IONIZING FLUXES DURING FIRST STAR EVOLUTION 



NLTE model atmospheres of first stars 



TABLE B2 
The same as Table IB1I but for 50Mq . 



age 






M* 






L 




l°gff 


Q(Hi) 


Q(Hc 


i) 


Q(Heii) 


r 

[years 




rn 1 

[Rq] 


[ M q] 






r r 1 


r t^i 

[K] 












[photons • s 1 













2.82 


50 


3 


597 


■ 10° 


84140 


5 


236 


2 


01 


10 48 


1.16 • 


10 48 


3 


80 


10 47 


9.203 


10 


2.99 


50 


3 


664 


■ 10 5 


82035 


5 


183 


2 


06 


10 48 


1.17 ■ 


10 48 


3 


79 


10 47 


2.545 


10 5 


3.08 


50 


3 


741 


■ 10 D 


81280 


5 


159 


2 


09 


10 48 


1.18 ■ 


10 48 


3 


82 


10 47 


4.606 


10 5 


3.17 


50 


3 


846 


• 10° 


80720 


5 


135 


2 


18 


10 48 


1.21 ■ 


10 48 


3 


93 


10 47 


6.683 


10 s 


3.27 


50 


3 


954 


• 10 5 


79980 


5 


107 


2 


21 


10 48 


1.23 ■ 


10 48 


3 


99 


10 47 


8.760 


10 D 


3.38 


50 


4 


064 


■ 10 5 


79250 


5 


079 


2 


27 


10 48 


1.26 • 


10 48 


4 


06 


10 47 


1.063 


10 6 


3.47 


50 


4 


178 


■ 10 5 


78700 


5 


055 


2 


34 


10 48 


1.29 • 


10 4S 


4 


00 


10 47 


1.432 


10 6 


3.71 


50 


4 


416 


■ 10° 


77270 


4 


999 


2 


48 


10 48 


1.34 • 


10 4S 


1 


07 


10 47 


2.024 


10 6 


4.18 


50 


4 


842 


• 10 5 


74470 


4 


895 


2 


72 


10 48 


1.44 • 


10 4S 


4 


32 


10 47 


2.359 


10 6 


4.56 


50 


5 


129 


• 10 5 


72280 


4 


818 


2 


87 


10 48 


1.47 ■ 


10 4S 


4 


23 


10 47 


2.642 


10 6 


4.99 


50 


5 


383 


■ 10 D 


69980 


4 


74 1 


3 


00 


10 48 


1.50 ■ 


10 48 


4 


21 


10 47 


2.957 


10 6 


5.63 


50 


5 


702 


■ 10 D 


66830 


4 


636 


3 


16 


10 48 


1.53 ■ 


10 48 


1 


18 


10 47 


3.289 


10 6 


6.68 


50 


6 


095 


■ 10 D 


62370 


4 


487 


3 


33 


10 48 


1.51 ■ 


10 48 


3 


89 


10 47 


3.584 


10 6 


8.20 


50 


6 


471 


• 10 5 


57150 


4 


309 


3 


11 


10 48 


1.43 • 


10 48 


3 


56 


10 47 


3.613 


10 6 


8.42 


50 


6 


516 


• 10 5 


56490 


4 


286 


3 


45 


10 48 


1.42 • 


10 4S 


3 


48 


10 47 


3.640 


10 6 


8.65 


50 


6 


561 


■ lO 5 


55850 


4 


263 


3 


46 


10 48 


1.40 • 


10 48 


3 


39 


10 47 


3.669 


10 6 


8.87 


50 


6 


592 


■ 10° 


55210 


4 


24 1 


3 


46 


10 48 


1.39 • 


10 48 


3 


33 


10 47 


3.866 


10 6 


11.1 


50 


6 


949 


■ 10 D 


50000 


i 


046 


3 


18 


10 48 


1.20 ■ 


10 4S 


2 


24 


10 47 


3.946 


10 6 


12.4 


50 


7 


046 


■ 10 D 


47530 


3 


952 


3 


11 


10 48 


1.09 ■ 


10 48 


1 


95 


10 47 


4.021 


10 6 


13.6 


50 


7 


228 


■ 10 D 


45600 


3 


869 


3 


39 


10 48 


1.02 ■ 


10 48 


1 


74 


10 47 


4.031 


10 6 


13.4 


50 


7 


311 


• 10 s 


46020 


3 


880 


3 


45 


10 48 


1.05 ■ 


10 48 


1 


82 


10 47 


4.032 


10 6 


13.2 


50 


7 


362 


• 10 5 


46560 


3 


897 


3 


51 


10 48 


1.09 • 


10 48 


1 


90 


10 47 


4.039 


10 6 


12.6 


50 


7 


621 


■ 10° 


48080 


3 


938 


3 


74 


10 48 


1.21 • 


10 48 


2 


18 


10 47 


4.047 


10 6 


12.8 


50 


7 


656 


■ 10 5 


47640 


3 


920 


3 


73 


10 48 


1.19 • 


10 48 


2 


13 


10 47 


4.064 


10 6 


13.1 


50 


7 


673 


■ 10 s 


47100 


3 


899 


3 


71 


10 48 


1.17 ■ 


10 48 


2 


06 


10 47 


4.075 


10 6 


13.5 


50 


7 


691 


■ 10 D 


46560 


3 


878 


3 


68 


10 48 


1.14 • 


10 48 


1 


99 


10 47 


4.084 


10 6 


13.8 


50 


7 


709 


■ 10° 


46020 


3 


857 


3 


66 


10 48 


1.11 ■ 


10 48 


1 


92 


10 47 


4.092 


10 6 


14.1 


50 


7 


727 


■ 10 D 


45500 


3 


836 


3 


63 


10 48 


1.09 • 


i n48 

10 


1 


86 


10 47 


4.099 


10 6 


14.5 


50 


7 


727 


• 10 5 


44980 


3 


816 


3 


60 


i r>48 

10 


1.06 • 


1Q 48 


1 


78 


10 


4.106 


10 


14.8 


50 


7 


745 


■ 10 


A A KdCX 

44060 


3 


799 


3 


58 


10 48 


1.04 • 


10 48 


1 


73 


10 47 


4.112 


10 6 


15.1 


50 


7 


762 


■ 10 5 


44050 


3 


778 


3 


55 


10 48 


1.01 • 


10 48 


1 


67 


10 47 


4.118 


10 6 


15.5 


50 


7 


780 


■ 10 5 


43550 


3 


757 


3 


52 


10 48 


9.83 • 


10 47 


1 


61 


10 47 


4.150 


10 6 


17.8 


50 


7 


852 


■ 10 5 


40740 


3 


637 


3 


31 


1Q 48 


8.25 ■ 


10 47 


1 


26 


10 47 


4.159 


10 6 


18.6 


50 


7 


870 


■ 10 5 


39900 


3 


600 


3 


23 


10 48 


7.77 ■ 


10 47 


1 


16 


10 47 


4.204 


10 6 


23.6 


50 


7 


980 


■ 10 5 


35480 


3 


390 


2 


73 


10 48 


4.33 ■ 


10 47 


3 


37 


10 46 


4.220 


10 6 


26.1 


50 


8 


017 


■ 10 5 


33810 


3 


304 


2 


48 


10 48 


3.14 ■ 


10 47 


2 


11 


10 46 


4.256 


10 6 


33.7 


50 


8 


139 


■ 10 5 


29850 


3 


081 


l 


75 


10 48 


9.08 • 


10 46 


3 


89 


10 45 


4.328 


10 6 


72.0 


50 


8 


355 


■ 10 5 


20560 


2 


422 


6 


34 


10 45 


2.04 • 


10 40 


3 


49 


10 35 


4.366 


10 6 


303. 


50 


8 


502 


■ 10 5 


10070 


I 


174 


I 


29 


10 41 


5.33 • 


10 35 


2 


07 


10 26 



Kubat 



TABLE B3 
The same as Table [Bl] but for 20Mq. 



a 


ge 












L 






Q(Hl) 


Q(Hei) 


Q(Heii) 


[years 




rn 1 

[Rq] 


[Af©] 






r r 1 


r t^i 

[K] 








[photons • s — 1 













1.65 


20 


4 


456 


■ 10 4 


65310 


5.305 


2.40 


10 47 


1.27 


-i ^A7 

10 47 


3.24 


10 46 


2 


720 


10 


1.73 


20 


4 


613 


■ 10 4 


64120 


5.258 


2.46 


10 47 


1.29 


10 47 


3.24 


10 46 


1 


581 


10 6 


1.90 


20 


5 


093 


■ 10 4 


62800 


5.179 


2.69 


10 47 


1.39 


10 47 


3.44 


10 46 


3 


164 


10 6 


2.15 


20 


5 


834 


■ 10 4 


61090 


5.072 


3.05 


10 47 


1.53 


10 47 


3.71 


10 46 


4 


726 


10 6 


2.52 


20 


6 


776 


• 10 4 


58610 


4.935 


3.52 


10 47 


1.70 


10 47 


3.91 


10 46 


6 


497 


10 6 


3.31 


20 


8 


317 


■ 10 4 


53830 


4.698 


4.13 


10 47 


1.80 


10 47 


3.49 


10 46 


7 


004 


10 6 


3.73 


20 


8 


913 


■ 10 4 


51640 


4.596 


4.30 


10 47 


1.79 


10 47 


3.32 


10 46 


7 


306 


10 6 


4.10 


20 


9 


480 


■ 10 4 


50000 


4.513 


4.45 


10 47 


1.79 


i nAI 

10 47 


3.23 


10 46 


7 


630 


10 6 


4.54 


20 


9 


772 


• 10 4 


47860 


4.424 


4.40 


10 47 


1.68 


i nAT 

10 47 


2.92 


10 46 


8 


012 


10 6 


5.34 


20 


1 


042 


• 10 5 


44870 


4.284 


4.33 


10 47 


1.52 


10 47 


2.49 


10 46 


8 


078 


10 6 


5.50 


20 


1 


057 


■ 10 D 


44360 


4.258 


4.32 


10 47 


1.49 


10 47 


2.42 


10 46 


8 


151 


10 6 


5.65 


20 


1 


074 


■ 10 D 


43950 


4.235 


4.34 


10 47 


1.48 


10 47 


2.38 


10 46 


8 


202 


10 6 


5.60 


20 


1 


096 


■ 10 D 


44360 


4.242 


4.50 


10 47 


1.55 


10 47 


2.51 


10 46 


8 


212 


10 6 


5.46 


20 


1 


132 


■ 10 D 


45290 


4.264 


4.78 


10 47 


1.69 


10 47 


2.78 


10 46 


8 


216 


10 6 


5.28 


20 


1 


191 


• 10 5 


46660 


4.294 


5.23 


10 47 


1.92 


10 47 


3.22 


10 46 


8 


235 


10 6 


5.09 


20 


1 


216 


■ 10 5 


47750 


4.325 


5.49 


10 47 


2.06 


10 47 


3.53 


10 46 


8 


267 


10 6 


4.94 


20 


1 


202 


• 10 5 


48300 


4.350 


5.49 


10 47 


2.09 


1 n47 
10 ' 


3.62 


10 46 


8 


405 


10 6 


5.36 


20 


1 


225 


■ 10 D 


46770 


4.286 


5.40 


10 


1.98 


10 47 


3.33 


10 


8 


506 


1U 


5.85 


20 


1 


250 


■ 10 


A A C TC\ 

448 tl) 


4.205 


5.25 


10 47 


1.83 


10 47 


2.95 


10 46 


8 


607 


10 6 


6.57 


20 


1 


279 


■ 10 5 


42560 


4.103 


4.99 


10 47 


1.61 


10 47 


2.49 


10 46 


8 


697 


10 6 


7.44 


20 


1 


312 


■ 10 5 


40270 


3.996 


4.65 


10 47 


1.38 


10 47 


2.03 


10 46 


8 


744 


10 6 


8.00 


20 


1 


334 


■ 10 5 


38990 


3.933 


4.47 


10 47 


1.27 


10 47 


1.82 


10 46 


8 


774 


10 6 


8.58 


20 


1 


349 


■ 10 5 


37760 


3.872 


4.20 


10 47 


1.13 


10 47 


1.58 


10 46 


8 


786 


10 6 


8.97 


20 


1 


358 


■ 10 5 


36980 


3.833 


4.02 


10 47 


1.04 


10 47 


1.43 


10 46 


8 


790 


10 6 


9.18 


20 


1 


358 


■ 10 5 


36560 


3.813 


3.87 


10 47 


5.84 


10 46 


5.18 


10 45 


8 


795 


10 6 


10.5 


20 


1 


312 


■ 10 5 


33960 


3.700 


2.81 


10 47 


1.97 


10 46 


1.25 


10 45 


8 


800 


10 6 


13.8 


20 


1 


318 


■ 10 5 


29580 


3.458 


1.11 


10 47 


1.68 


10 44 


1.58 


10 41 


8 


805 


10 6 


19.9 


20 


1 


297 


■ 10 5 


24550 


3.141 


5.78 


10 45 


7.50 


10 41 


1.91 


10 38 



NLTE model atmospheres of first stars 



TABLE B4 
The same as Table IB1I but for IOMq . 



age 












L 


T c ff 


l°gff 


Q(Hi) 


Q(Hc 


i) 


Q(Heii) 


r 

[years 




[-R©] 


[Af©] 






r r 1 

[£© 


[K] 










[photons • s 1 













1.37 


10 


6 


653 


■ 10 J 


44460 


5 


162 


2.48 


10 46 


8 


36 • 


10 45 


1 


39 


10 45 


4.445 


10 


1.38 


10 


6 


887 


• 10 3 


44770 


5 


159 


2.61 


10 46 


8 


94 • 


10 45 


1 


50 


10 45 


3.584 


10 6 


1.39 


10 


7 


852 


■ 10 d 


46020 


5 


150 


3.12 


10 46 


1 


13 ■ 


10 46 


1 


94 


10 45 


5.616 


10 6 


1.41 


10 


8 


690 


■ 10 J 


46990 


5 


142 


3.57 


10 46 


1 


34 ■ 


10 46 


2 


34 


10 45 


8.502 


10 6 


1.47 


10 


9 


954 




47530 


5 


103 


4.15 


10 46 


1 


59 ■ 


10 46 


2 


82 


10 45 


1.068 


-i r\7 

10' 


1.62 


10 


1 


081 


■ 10 4 


46240 


5 


019 


4.42 


10 46 


1 


63 • 


10 46 


2 


83 


10 45 


1.296 


10' 


1.83 


10 


1 


197 


■ 10 4 


44560 


1 


911 


4.60 


10 46 


1 


62 • 


10 46 


2 


76 


10 45 


1.493 


-i n7 

10' 


2.13 


10 


1 


327 


■ 10 4 


42460 


4 


782 


4.73 


10 46 


1 


56 • 


10 46 


2 


57 


., „ ,1c. 

10 45 


1.658 


i n7 

10 


2.52 


10 


1 


469 


■ 10 4 


39990 


4 


634 


4.51 


10 46 


1 


34 • 


10 46 


2 


14 


10 45 


1.725 


10' 


2.76 


10 


1 


541 


■ 10 4 


38720 


4 


557 


4.20 


10 46 


1 


13 ■ 


10 46 


1 


79 


10 45 


1.746 


-i r\7 

10' 


2.84 


10 


1 


567 


■ 10 4 


38280 


4 


530 


4.11 


10 46 


1 


09 ■ 


10 46 


1 


71 


10 45 


1.769 


-i r\7 

10' 


2.92 


10 


1 


596 


■ 10 4 


37930 


4 


506 


4.06 


10 46 


1 


06 ■ 


10 46 


1 


67 


10 45 


1.778 


-i r\7 

10' 


2.95 


10 


1 


614 


■ 10 4 


37840 


4 


497 


4.08 


10 46 


1 


05 ■ 


10 46 


1 


65 


10 45 


1.791 


i r\7 
10' 


2.93 


10 


1 


652 


■ 10 4 


38190 


4 


503 


4.32 


10 46 


1 


14 • 


10 46 


1 


80 


10 45 


1.794 


i r\7 

10' 


2.89 


10 


1 


679 


■ 10 4 


38640 


4 


516 


4.57 


10 46 


1 


24 • 


10 46 


1 


96 


10 45 


1.795 


10' 


2.87 


10 


1 


714 


■ 10 4 


38990 


4 


523 


4.81 


10 46 


1 


34 • 


10 46 


2 


12 


10 45 


1.795 


f r\7 

10' 


2.84 


10 


1 


750 


■ 10 4 


39350 


4 


530 


5.06 


10 46 


1 


44 • 


10 46 


2 


28 


10 45 


1.796 


i r\7 

10' 


2.82 


10 


1 


791 


■ 10 4 


39720 


4 


536 


5.33 


10 46 


1 


56 • 


10 46 


2 


51 


10 45 


1.796 


10' 


2.80 


10 


1 


832 


■ 10 4 


40090 


4 


542 


5.60 


10 46 


1 


67 ■ 


10 46 


2 


68 


10 45 


1.796 


i / i 7 

10' 


2.79 


10 


1 


879 


■ 10 4 


40460 


4 


547 


5.89 


10 46 


1 


79 ■ 


10 46 


2 


88 


10 45 


1.804 


i r\7 

10' 


2.74 


10 


1 


950 


• 10 4 


41210 


4 


563 


6.41 


10 46 


2 


02 • 


10 46 


3 


27 


10 45 


1.810 


-i r\7 
10' 


2.68 


10 


1 


932 


• 10 4 


41590 


4 


583 


6.49 


10 46 


2 


07 • 


10 46 


3 


36 


10 45 


1.825 


-i n7 
10' 


2.61 


10 


1 


928 


■ 10 4 


42070 


4 


604 


6.64 


10 46 


2 


16 • 


10 46 


3 


51 


10 45 


1.829 


10' 


2.61 


10 


1 


928 


■ 10 4 


42070 


4 


604 


6.64 


10 46 


2 


16 • 


10 46 


3 


51 


10 45 


1.862 


i n7 
10' 


2.68 


10 


1 


963 


■ 10 4 


41690 


4 


580 


6.63 


10 46 


2 


13 ■ 


1 n46 

10 


3 


16 


10 45 


1.878 


10' 


2.75 


10 


1 


986 


■ 10 4 


41300 


4 


559 


6.58 


10 ™ 


2 


07 • 


10 46 


3 


35 


i n45 

10 


1.888 


10 


2.83 


10 


2 


009 


i n4 
• 10 


a no on 
4083U 


4 


534 


6.47 


10 46 


2 


00 ■ 


10 46 


3 


21 


10 45 


1.899 


10 7 


2.91 


10 


2 


042 


■ 10 4 


40460 


4 


511 


6.56 


10 46 


2 


01 ■ 


10 46 


3 


20 


10 45 


1.909 


10 7 


2.98 


10 


2 


075 


■ 10 4 


40090 


4 


488 


6.55 


1Q 46 


1 


99 • 


1Q 46 


3 


17 


10 45 


1.919 


10 7 


3.06 


10 


2 


109 


■ 10 4 


39720 


4 


465 


6.50 


1Q 46 


1 


94 • 


10 46 


3 


07 


1Q 45 


1.927 


10 7 


3.16 


10 


2 


143 


■ 10 4 


39260 


4 


438 


6.40 


10 46 


1 


87 • 


10 46 


2 


93 


1Q 45 


1.938 


10 7 


3.93 


10 


2 


000 


■ 10 4 


34590 


4 


248 


3.25 


10 46 


1 


43 • 


10 45 


9 


46 


10 43 


1.940 


10 7 


5.25 


10 


2 


104 


■ 10 4 


30340 


3 


998 


8.90 


10 45 


1 


08 ■ 


10 43 


8 


21 


10 39 


1.943 


10 7 


11.4 


10 


1 


946 


■ 10 4 


20180 


3 


324 


8.52 


10 43 


1 


04 • 


10 39 


3 


32 


10 34 


1.943 


10 7 


12.3 


10 


1 


928 


■ 10 4 


19360 


3 


256 


7.27 


10 43 


9 


12 • 


10 38 


6 


10 


10 33 


1.943 


10 7 


15.9 


10 


1 


897 


■ 10 4 


16980 


3 


035 


2.04 


10 43 


1 


75 ■ 


10 40 


1 


61 


10 34 



